Einfaches Triebwerksmodell zur Simulation nichtlinearer 
dynamischer Schubantworten 

Günther SCHÄNZER und Peter KRAUSPE, TU Braunschweig *) 

Übersicht: Es wird eine einfache mathematische Triebwerksbeschreibung gegeben, mit der gemessene Übergangsfunktionen des 
Triebwerksschubes sowohl im quasistationären Betrieb als auch im Beschleunigung*- bzw. Verzögerungsbereich nachgebildet 
werden können. Das Modell ist aufgrund seiner Einfachheit hei gleichzeitig hoher Nachbildungsgüte besonders zur Darstellung 
von TL-Triebwerken in digitalen Simulationsprogrammen geeignet. 

Simple mathematical model for simulating nonlinear jet engine thrust response 

Summary: A model for the Performance ofjet engines is introduced which allows to shnulate measured thrust-to-time curves in 
steady-state Operation as well as in accelerationi 'deceleration scbedale. Due to its simplicity along with its high fidelity it is very 
well suited for the representation of jet engines in digital Simulation pragrams. 
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1, Einleitung 

Die Darstellung von gemessenen Schub-Zeit- Vet laufen von 
Strahltriebwerken in digitalen Simulationsprogrammen, an 
deren Genauigkeit besonders dann hohe Anforderungen 
gestellt werden, wenn der Triebwerksschub die primäre 
Stellgröße eines Regelungsprozesses ist, gestaltet sich oft 
problematisch. Zum einen lassen sich die Übergangs- 
funktionen (BILD 1) aufgrund der großen Verhältnisse von 
Verzugszeit und Ausgleichszeit nur durch lineare Regel- 
strecken sehr hoher Ordnung (20 und mehr) annähernd 
darstellen. Zum anderen verändert sich das Ubergangs- 
verhalten des Schubes nichtlinear in Abhängigkeit von der 
Triebwerksdrehzahl (BILD 2). Schließlich ist das Zeit- 
verhalten durch trieb werks interne Begrenzungen des zu- 
lässigen Luftmassenstroms von der Größe des Stell- 
größensprunges abhängig, wobei zwischen dem Betrieb auf 
der sogenannten steady-state operating line (quasistatio- 
närer Betrieb) und der sogenannten acceleration/decelera- 
tion schedule (Beschleunigung bzw. Verzögerung) unter- 
schieden werden muß. Die Anpassung einer linearen Regel- 
strecke hoher Ordnung an die genannten Bedingungen ist 
sehr aufwendig [l]. Die Beschreibung der Triebwerks- 
dynamikinder Phasenebene [2] bietet hier die Möglich- 
keit einer einfachen und übersichtlichen Darstellung der 
linearen und nichtlinearen Anteile. Die verwendeten 
Modellparameter sind abhängig vom Triebwerks- und 
Flugzustand, wie z. B. Druckverhältnisse, Luftmassen- 
durchsatz, Einlauftemperatur, Machzahl u. ä. Aus Gründen 
der Einfachheit witd bewußt darauf verzichtet, die physi- 
kalischen Zusammenhänge zwischen den Triebwerks- 
parametern mit in das Modell einzubeziehen, wie es sonst 
bei der Simulation des gesamten Triebwerks üblich ist [3]. 
Die Modellparameter des vorgestellten Modells werden 
grundsätzlich aus Schubmessungen {Bodenstandsversuch 
oder Simulation) gewonnen, die für den jeweiligen Flug- 
zustand gültig sind. 
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BILD 1 : Obergangsverhalten verschiedener Strahltriebwerke 



JOB 

% 
F 80 

SB 
L0 

28 

2 U S 8 s W 

t 

BILD 2: Sprungantworten eines Strahltriebwerks (Messungen) 
für a) große und b) kleine Schubänderungen 
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BILD 3 : Phasenporträt des Schubverlaufs eines Strahltriebwerks 



3. Grundlagen und physikalische Zusammenhänge 

Das dynamische nichtlineare Triebwerksverhalten, wie 
es z.B. in BILD 2 wiedergegeben ist, läßt sich zumindest 
für die untersuchten Mehrwellentriebwerke JT3D.JT9D, 
CF6-50C in der Phasenebene F,F mit guter Genauigkeit 
durch Geraden approximieren (BILD 3). Die gesuchte 
Funktion F (t) kann, da jedem Wert F ein Wert F zugeordnet 
werden kann, aus einer einfachen Integration von F ge- 
wonnen werden: 



(1) 



F{t)=F + !P(F)dt. 
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BILD 4: Vergleich gemessener und simulierter Schubverläufe 
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In BILD 4 ist ein Vergleich der auf diese Weise berechneten 
Schub-Zeit- Verläufe mit den zugrunde liegenden Messungen 
aufgetragen. Durch die zusätzliche Angabe der Drehzahl Ni 
des Hochdruckverdichters treten die typischen dynamischen 
Eigenschaften eines Strahltriebwerks deutlicher hervor 
(BILD 5). Unterteilt man die Phasenkurven in BILD 3 und 
die zugeordnete Sprungantwort auf eine Gasdrossel- 
änderung in BILD 5 in Teilbereiche, so läßt sich für jeden 
Einzelbereich eine anschauliche physikalische Deutung 
finden: 

(I) Anlaufbereich 

Hochlaufen der Triebwerksdrehzahl Ni aus dem Leerlauf 
mit nahezu konstanter, relativ kleiner Drehzahländerung 
Ni, [, Es ist zu beachten, daß bei Turbinen-Luftstrahl-Trieb- 
werken ein starker nichtlinearer Zusammenhang zwischen 
dem Schub und der Drehzahl besteht. So ist z. B, dem Leer- 
laufschub F L eine Leerlaufdrehzahl von fast 60% N2 vmax 
zuzuordnen; innerhalb der restlichen 40% des zulässigen 
Drehzahlbereichs liegt somit das gesamte Schubspektrum 
von Fl bis fma.v Im Bereich 1 ergibt sich trotz des nicht- 
linearen Verhaltens der Drehzahl ein nahezu konstanter 
Wert F. 

(II) Übergangsbereich auf maximale Dreh- 
beschleunigung N2, raJ ^ 

[n der Phasenebene ist die Steigung F/F der Phasenkurve 
konstant. 

(III) Bereich konstanter Schub ander ung F wix 
mit fast linearer Abhängigkeit von 
Schub und Drehzahl 

Die Beschleunigung eines Mehrwellentriebwerks wird in 
diesem Bereich durch die Pumpgrenze des Hochdruck- 
verdichters begrenzt. Bei Mehrwellentriebwerken w-ird 
der sichere Abstand zur Pumpgrenze durch Trieb werks- 
regler gewährleistet. Damit wird das dynamische Verhalten 
und besonders das maximale Beschleunigungsvermögen 
im Bereich III durch das Zusammenwirken von Triebwerk 
und Regler bestimmt. 

(IV) Übergangs be reich auf den stationären 
Endwert 

Das Zeitverhalten entspricht dem Verhalten eines linearen 
Systems erster Ordnung mit der charakteristischen Zeit- 
konstante Tj, (BILD 6): 



(2) F{t) = -^-[F*a(t)-FM]= -—AF(t). 



Die Steigung der Phasenkurve in der Phasenebene ergibt 

sich daraus zu 















i 


-Wh F m „ 








h 














*N 












-SO 

-so ED . 


'm 






N z IFi | 












-CO 

, II 


1, / 














-26 J 




N l 






1" " "" 




f 










fis: -ta __- 


100% SS 
■■:-„ 


SS LS 


2S 




X 2 
-2 \ 




L 


s 




! 5 


■V, 














■%* 




-i 


1 




1/ 


N,lt) 


f | JT3B 
















1 


- 


















tl 













BILD 5; Abhängigkeiten von Schub und Drehzahl 
eines .Mehrwellentriebwerks 
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Ein Vergleich mit BILD 5 zeigt, daß bei dem hier unter- 
suchten Triebwerk näherungs weise 



(5) 



T 1 =T 1 = T, 



gesetzt werden kann. Diese Zeitkonstanten charakteri- 
sieren das Zeitverhalten des Triebwerks bei kleinen Soll- 
wertsprüngen {„Kleinzeitverhalten") und liegen im Bereich 
von 0,25 bis 1 sec. Es ist jedoch anzumerken, daß die Modell- 
parameter E max , T„ T] und Ti jeweils nur für einen bestimm - 
ren Umgebungs- und Flugzustand gelten. Die hier ange- 
gebenen Modellparamerer beziehen sich auf Bodenstands- 
versuche unter Normalatmosphärenbedingungen. Wie aus 
Flugmessungen bekannt ist, ergeben sich in größerer Höhe 
infolge der abnehmenden Luftdichte und des demzufolge 
kleineren Luftmassenstroms größere Werte für T] und T 2 . 

Mit Hilfe des Kleinzeirverhaltens lassen sich auch jene 
quasilinearen Zustandsänderungen beschreiben, die durch 
relativ kleine Stellgrößensprünge gekennzeichnet sind, 

denen das System ohne große Verzögerungen folgen kann. 
Die charakteristische Zeitkonstante T„ des linearen 
Systems kann durch Tangentenbildung direkt dem ge- 
messenen Verlauf der Übergangsfunktion der Strecke 
entnommen werden, und zwar, wie BILD 2 zeigt, aus der 



(3) 



( P\ 1 

\AF/iv T, 



Die charakteristische Zeitkonstante des Bereichs IV sei T 2 . 
Vergleicht man die Übergangsbereiche II und IV, so läßt sich 
eine formale Übereinstimmung finden: Beide Abschnitte 
werden in der Phasenebene durch eine Gerade mit der 
Steigung F/F wiedergegeben. Im Bereich II sei analog zu 
Gl. (3) die charakteristische Zeitkonstante T t definiert: 





BILD 6: Spmnganrwon und Phasenporträt 
eines linearen Systems erster Ordnung 
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Sprungantwort sowohl auf a) große als auch b) kleine 
Sprünge; 



(6) 



-T b =T v . 



Wesentlich ist, daß die Strecke praktisch alle möglichen 
Zwischenwerte von F annehmen kann, die unterhalb des 
maximalen Wertes f max der Beschleunigungslinie liegen. 
In der Phasenebene wird das quasilineare Verhalten des 
Triebwerks demnach durch die Menge aller Punkte (F,F) 
innerhalb der Phasenkurve wiedergegeben. 




BILD 7: Blockschaltbild des Triebwerkmodells 



3. Modellbeschreibung 

Die beschriebenen Zusammenhänge sind im Blockschalt- 
bild in BILD 7 dargestellt. Nach Gl. (2) wird in Abhängigkeit 
von der Differenz zwischen Soll- und Istwert des Schubes F 
der Anfangswert F eines linearen Systems mit der Zeit- 
konstanten TV ermittelt. Liegt der dadurch bestimmte Punkt 
(F,F) in der Phasenebene innerhalb der durch die Phasen- 
kurve begrenzten Fläche (BILD 8 a), so handelt es sich um 
eine quasilineare Zustandsänderung, und der gesuchte 
Wert F ist gleich dem Ergebnis von Gl. (2); liegt der Punkt 
außerhalb der Phasengrenzen (BILD 8 b), so wird für F die 
maximal mögliche Schubänderung entsprechend dem zuge- 
hörigen Ordinatenwert der Phasenkurve eingesetzt. In 
BILD 7 wird dieser Vorgang durch die Funktion des Be- 
grenzers B realisiert, dessen Kennlinie die Phasenkurve 
des Triebwerks ist. Eine anschließende Integration von F 
nach GL (1) liefert den gewünschten Wert F{t). Entspre- 
chendes gilt für den Betrieb auf der Verzögerungslinie 
bzw. im Bereich des quasistationären Verzögerns des 
Triebwerks. 

Die gute Übereinstimmung der gemessenen und simulierten 
Schubverläufe, verbunden mit dem geringen Programmier- 
aufwand, lassen das vorliegende Triebwerksmodell als 
sehr brauchbar erscheinen. Anwendungen mit detaillierten 

Zahlenrechnungen sind in [4] enthalten. 
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BILD 8 : Beispiele von Schubkommandos in der Phasenebene 

a) 75— »80% (quasilinear), 
b) 40—» 50 % (begrenzt beschleunigt) 
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